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Abstract. We address the quantum-classical correspondence for chaotic systems 
with a crossover between symmetry classes. We consider the energy level statistics 
of a classically chaotic system in a weak magnetic field. The generating function of 
spectral correlations is calculated by using the semiclassical periodic-orbit theory. An 
explicit calculation up to the second order, including the non-oscillatory and oscillatory 
terms, agrees with the prediction of random matrix theory. Formal expressions of the 
higher order terms are also presented. The nonlinear sigma (NLS) model of random 
matrix theory, in the variant of the Bosonic replica trick, is also analyzed for the 
crossover between the Gaussian orthogonal ensemble and Gaussian unitary ensemble. 
The diagrammatic expansion of the NLS model is interpreted in terms of the periodic 
orbit theory. 
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1. Introduction 

The quantum-classical correspondence in a chaotic system is a longstanding problem. 
In a quantum system, chaos does not exist in the sense of the "sensitivity to the initial 
conditions", while in a classical system, it exists and a characterization is given by a 
positive Lyapunov exponent. A signature of chaos in a quantum system can be seen in 
the energy level statistics, which shares the same behavior with the eigenvalue statistics 
of random matrices [1]. According to a well-known conjecture by Bohigas, Giannoni, 
and Schmit |2|, the energy level statistics in classically chaotic systems with and without 
time-reversal symmetry are universally described by the Gaussian orthogonal ensemble 
(GOE) and Gaussian unitary ensemble (GUE), respectively. 

Recently a justification of this conjecture has been given for the spectral form 
factor, in terms of the semiclassical periodic-orbit theory p]. It is based on Berry's 
diagonal approximation |3] for the first order term and the periodic-orbit pairs with an 
encounter [5] for the second order term. Higher order terms are also calculated by using 
the periodic pairs with many encounters. If a very weak magnetic field is applied to a 
classically chaotic system, a crossover between the GOE and GUE universality classes 
is realized. The level statistics in this domain is described by Pandey and Mehta's two 
matrix model [6]. The periodic-orbit theory is also applicable to explain the behavior 
in this GOE-GUE crossover domain [71 El El [IDl El- In this approach, a stochastic 
behavior of the periodic orbits plays a crucial role. Similar ideas also give evaluations 
of GUE-GUE[ni[12], G0E-G0E[12J, GSE-GSE[I3] and G0E-GSE|I3] form factors. 

Moreover a recent progress on this topic has made it possible to reproduce the 
oscillatory terms in the expansion of the GOE and GUE spectral correlations [14J. In 
this new scheme, a periodic orbit theory is constructed for the generating function which 
yields the spectral correlation. The generating function is defined as 

/ det{E^-H) det {E-^-H) \ 
Z{e^, es, ec, en) - ^ _ _ j . (1) 

where H is the Hamiltonian. Here E^^j and E]^ are given by Ej^^ = E + and 

E]^ ]j = E + with the local average of the level density p, where the offsets from the 
center energy E are taken with an imaginary part as = e ± Z7 (7 > 0) to ensure the 
convergence. The generating function yields the spectral correlation function C(e) 
as 

d^Z 1 „ / 1 \ „ / 1 



deAdes II (2vrp)2^' [e + ^. + - H ) ^\e - - - H ^ 

= lc{e) + \, (2) 

where the symbol || denotes the identification e\ = e'^ = e'*",e^ = ~ ^ crucial 
idea of the new scheme is counting all the diagrams in terms of the pseudo orbits, which 
correspond to the terms appearing in a field theoretical treatment of random matrices. 
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In this paper, we extend this new scheme to the case with a weak magnetic field, 
which describes the GOE-GUE crossover. We combine the ideas of Refs.[Tl] and |10j . 
and compare the results of periodic-orbit and random matrix theories. The random 
matrix theory prediction for the generating function is known to be |15j . 

Z{eA,eB,ec,eD) = Z'^'^ + Z^'\ (3) 



where 



(e^ — eB)(ec — cd) 
= dp p-' e^b-i)^AB/2-2sp2 y ' dq q e-^(7-i)ecz,/2+2V 

2 2 2 2 

+ 



iieAB + iSs) i{eAD + i8s) i(ecB + «8s) i(ecz) + «8s) 

1 

l6seAC(^BD 



{ecD + i8sy{eAD + i8s){ecB + i8s 
1 



_^ C5) 

{eAB + ^8sy{eAD + i8s){ecB + ^8s) 
with €ai3 = ea — G/s- The function Z^'^^ is obtained from the function Z^^^ by exchanging ec 
and eo- The parameter s controls the symmetry: the GOE (GUE) limit corresponds to 
s — J- (s — )■ oo). In this paper, we explicitly calculate the generating function up to the 
second order using the periodic-orbit theory, and confirm the agreement with the above 
prediction. Moreover, in order to see the correspondence of the higher order terms, we 
examine the nonlinear sigma (NLS) model in the crossover domain, and interpret the 
diagrammatic expansion of the NLS model in terms of periodic pairs. 

2. Semiclassical Expression of the Generating Function 
2.1. Magnetic Action 

We consider a bounded quantum system with / degrees of freedom whose corresponding 
classical dynamics is chaotic. We assume that the system has a gauge potential due to 
a magnetic field B. Let us denote the energy by E and each phase space point by a 2/ 
dimensional vector x = {q,p), where / dimensional vectors q and p specify the position 
and momentum, respectively. In the semiclassical limit /i — )> 0, the trace formula for the 
energy-level density implies 

Tr {E+ -H)-' ^ - mp - T,F,eH[^'^(^+)+©'^(^)l, (6) 

a 

where H is the Hamiltonian. The first term in the right hand side gives the local average 
of the level density p and the second term describes the fiuctuation around the average. 
The local average of the level density is equal to the number of Planck cells inside the 
energy shell 

P = (2^' 
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where Q{E) is the volume of the energy shell. The function Sa is the classical action 
(including the Maslov phase) of the periodic orbit a with a period = dSa/dE, Fa is 
the stability amplitude and 6a is the magnetic action. 

We consider the following two possible situations. One is the case that the system 
has an Aharonov-Bohm type gauge potential, where the gauge potential exists although 
the magnetic field is not directly applied to the system. Then the classical dynamics is 
time- reversal invariant, while the quantum counterpart is not. Another situation is the 
case that the magnetic field is directly applied to the particle, but it is sufficiently weak 
such that the cyclotron radius is much larger than the system size and thus the presence 
of the magnetic field does not significantly change Q{E). In this case the impact of 
the magnetic field on the orbits can be neglected. In both cases the magnetic field 
still influences the behavior of the system through its contribution to the phase in (6). 
Time-reversed classical dynamics must be taken into account to derive the form factor 
in the crossover domain; it even becomes nontrivial to derive the form factor in the GUE 
limit. 

The magnetic action 0a(-B) is a function of the magnetic field and is defined as 

dqa 



QaiB) = B / A{q) ■dq = B gait) dt, g^it) = A^q^ 



dt 



(8) 



where A{q) is the gauge potential which generates the unit magnetic field and qaif) 
describes a classical motion in the configuration space along the orbit a. Since the 
classical dynamics is chaotic, successive changes of the velocity can be regarded as 
independent events. Hence, if the time T elapsed on an orbit a is sufficiently large, the 
statistics of ga{t) can be regarded as that of the Gaussian white noise satisfying the 
correlation {{gait) gait'))) = 2D6it — t'). The Gaussian average ((■ ■ ■)) is computed as a 
functional integral 



(( F[ga] )) 



j Vga exp 








j Vga exp 







(9) 



It is known that this averaging yields the non-oscillatory terms of the spectral 
correlation [lOl [11] . 



2.2. Pseudo-Orbit Representation of the Generating Function 
Utilizing Eq.dH]), we can rewrite det(-E'+ — H) as 

det(E+ - if) oc exp |y dE' Tr ^ | 

oc exp -mAr(E+)-^Fae*^^"(^^)+®"(^» 



(10) 
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Expanding the exponential we obtain an equivalent expression in terms of the pseudo 
orbits [16]: 

det{E+ - if) oc e-*-^(^^) 1 1 _ ^ F„e^(^"+®") 

a,a' J 



-iiTN{E+) 



(£;+)+0a(b)} 



:iii 



where A is the index of a pseudo orbit, which consists of component orbits, and 
is the number of the component orbits. The functions Sa_ and Byi are the sums of the 
mechanical and magnetic actions of the component orbits, respectively, and F4 is the 
product of the stability amplitudes. The corresponding factor det(i?~ — H) is given by 
the complex conjugation of det (£'"'" — H) as 

det(E- -H)oc e-^(^-) ^F;(-l)^-^e-^^^^(^")+®^(^», (12) 

A 

where an asterisk means a complex conjugate. 

Berry and Keating [16] postulated that Eq. (fTT]) and Eq. (fT2|) become identical in 
the limit E^^E~ — )■ E and found a duality relation between the contribution of the 
pseudo-orbits with a duration larger than a half of the Heisenberg time Th = 27rhp and 
those with a duration shorter than Th/2. Thus the sum over the long orbits can be 
treated as the complex conjugate of the sum over the short orbits. This leads to the so 
called Riemann-Siegel lookalike formula 

det(E - if) oc e-'-^(^) J2 i^A(-l)^^e^^^^(^)+®-^(''» + c.c., (13) 

Ta<Th/2 

This formula has an analogy in the corresponding expression of Riemann's zeta function. 
Using Eqs. OlOl) and ffTTl) . we write the determinants of the inverse matrices as 



det{E^ 



H 



-1 



oc exp 



z7riV(i5+) + 5^i^e*l^''(^^)+®" 



(14) 



det(E" - Hy^ oc exp 



-l7,N{E-) + Y.^a 



e A 



{Sa{E-) + ea{B)} 



)+^a{B)} 



(15) 



Inserting these expressions directly into the definition of the generating function, we 
obtain a semiclassical expression : 



-iiTN(Ec) 



E ■ 

Tc<Th/2 
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Td<Th/2 

l{SA{E+)+eA} 



B 



(16) 



We used the Riemann-Siegel lookalike formula (fT3|) and its complex conjugate for the 
numerator and the expressions fll4p - (IT3]) for the denominator. Because the expression 
f fT3|) has two parts, Eg. (1161) in total has four parts, each of which has the factor 

^±i-K[N{E+eA/{2TTp)}-N{E-eB/{2-Kp)}\ ^±iTT[N {E+tA/ {2ivp)}+N {E-eg / {2Tvp)}\ _ HoWCVer Cach part 

with the latter factor has no contribution because it is expected to vanish after averaging 
over E due to a rapid oscillation. 

Eventually we obtain two parts for Z as 



(17) 



where the first term Z^^^ is written as 



zm 



E 



oii^\-^B-''C+^D)/2 



A,B,C,D 
(Tc,Tr,<TH/2) 



p^p*^p^p*^i^^l)Nc+ND^i{SA{E)-SB{E)+Sc{E)-SD{E))\ 



I 



Here the symbol ((...)) implies the average over the weight Eq.([9]). To derive Eq. ffTSj) . 
we used the expansion N{E ± ejlup) ~ N{E) ± e/27r,and S{E ± e/27rp) jh ~ S{E) jh ± 
TejTii. The term Z'^'^^ is obtained from Z^^^ by exchanging the variables ec and as 

) Cfi) CC) ^d) — (19) 

Utilizing these expressions, we write the second derivatives of Z, which gives the spectral 
correlation C(e), as 



d^z 



deAdes 



deAdes 



52^(2) 



deAdee 



deAdes 



+ 



deAdee 



(20) 



where x denotes the identification 



7 +0; CD 



e. 



In 



Ref. 113], the identifications with the symbols || and x were called "column- wise" and 
"crosswise" , respectively. 



3. Consecutive Approximation Orders in Terms of Periodic Orbits 

3.1. Diagonal Approximation 

Let us first consider the diagonal approximation, in which different orbits are assumed 
to be uncorrelated, except for the orbits obtained from each other by time reversal. We 
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express Eq. (IT5i) in terms of the periodic orbits, not using pseudo orbits. On account of 
Eq. ([TOl) . the generating function Z^^^ can be written in terms of the periodic orbits: 



AC 



_f'a * 
JBD 



(21) 



It can be rewritten as an averaged product over the pairs (a, a'^^) of the mutually 
time reversed orbits; considering that a time-reversed orbit has a different sign in the 
gauge-potential term, we obtain 



^j(e^-eB-ec+ei3)/2 




-p) \ 



(22) 



Now let us calculate the average. As the factors in the product are assumed to be 
uncorrelated, the averaged product of the factors can be replaced by the product of the 
averaged factors. The average over the energy E is replaced by an average 

over the classical action Sa- The terms containing the magnetic action are averaged over 
the Gaussian white noise (see Eq.(l9]) as 

((e^))=e-^^/^-, (24) 

where h = AB'^DTfj/h^- Considering that the coefficients flcy fsD ^'^^ exponentially 
small for the relevant long orbits, we expand each factor up to the second order and 
doubly average the result. Then we find 



(^exp 



.iSa fa 



a ea * 



AC\ 



7, 



BD 



eft " + e 



1 + /ac/V (^2 + ((enr)) + ((e^)) 

^exp [2/1^/V (l + e-^^»/^-)] 
It is now straightforward to obtain the expression 

7^^^ ~ pi(eA-«i3-«C+<;D)/2 
^diag ^ 



(25) 



X exp 



(26) 



where summation in the exponent is carried over orbits, not orbit pairs; for 
compensation, the factor 2 is dropped. 
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We can further proceed with the calculation by using Hannay and Ozorio de 
Almeida (HOdA) sum rule [H] 



FaW-) - I ^(O, (27) 



where the lower limit of the integration is the minimum period Tq; in the semiclassical 
limit it can be replaced by zero. Then we immediately obtain the result of the diagonal 
approximation: 

_ „j(€A-es-ec+eD)/2 

diag 

^ [ {ec - ei?)(eA -ep) } f (eg - + ib){eA - ep + ih) \ 
\{eA-eB){ec-ep)]\{eA-eB + ib){ec-ep + ib)]' 

3.2. Second Order Contribution 

Let us now go beyond the diagonal approximation. In order to calculate the off-diagonal 
contributions, we can utilize the method developed in [21 E]. In that method, the 
periodic orbits are divided into links and encounters. An encounter is a part of an 
orbit where different links meet together, and orbit stretches connect the links inside 
the encounter region. It is assumed that the contributions come from the periodic-orbit 
pairs with the same (or time-reversed) links which are differently connected inside the 
encounter regions. Due to the different connections in the encounter regions, there is a 
small but significant difference in actions between the original and partner orbits. The 
gauge potential acts on both links and encounters. The gauge potential constructively 
interferes when the original and partner orbits go in the opposite directions. On the 
other hand, when the directions are the same, a destructive interference takes place. 

The off-diagonal contribution is formally written by using the pseudo-orbit 
representation (fT8|) as 



A,B,C,D 

where YI'a b c d i^^ans a summation over the pseudo orbits forming links and encounters. 
AS" represents a sum of action differences in the encounter regions, and AB expresses 
the effect of interferences of the gauge potential. The overall sum Z^^^ is obtained by 
multiplying the two contributions 

Z^'' = ZZZS^. (30) 

Let us now calculate the leading (second order) term of the off-diagonal 
contribution. We first discuss the Sieber-Richer (SR) pairs, one of which is depicted 
in Figure 1. A SR pair has one encounter and two links, and the original and partner 
orbits go in the opposite directions in one of the links. In Figure 1, the solid curve is 
the original orbit and the dashed curve is the partner. A pseudo-orbit index "A" or "C" 
is assigned to the original orbit, and "B" or "D" is to the partner. There are in total 
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Figure 1. An example of the Sieber-Richter (SR) pairs, where "A" is assigned to the 
original orbit (the solid curve) and "B" is to the partner (the dashed curve). 
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Table 1. A schematic table for the 16 possibilities of the SR pairs. Here means the 
absence in the diagram. 



16 possibilities, as listed in Table 1. On the Poincare section within the encounter 
region, the phase space can be furnished with canonical coordinates u and s, which 
are the coordinates in the unstable and stable directions, respectively. Let us suppose 
that the original orbit pierces the Poincare section at the origin (0, 0) and at the point 
{u,s). To demand that the piercing points in the encounter are mutually close, we 
introduce a bound c and assume that |s| < c. Then we can estimate the duration 
tenc of the encounter and the action difference AS of the original and partner orbits 
as tenc = i 1^ ]^ and AS* = us, where A is the Lyapunov exponent. Now we need to 
estimate the number of encounters in one periodic orbit of a period T. This can be 
computed as 

fT-2U 



duds I dt J' (31) 

where Qr^ ^= 2tt\th ) ergodic return probability per unit action. The 

combinatorial factor 2^ implies that one of the two encounter stretches is chosen as 
the first stretch and one of the two directions is chosen as the positive direction. The 
factor T "^^""^^ dt indicates that one of the two piercings occurs in the time interval 
[0, T] and that the time t elapsed on one of the links lies in [0, T— 2tenc]- The contribution 
from the gauge potential has to be taken into account on one of the links, where the 
original and partner orbits go in the opposite directions. In the SR pair, a constructive 
interference of the gauge potential does not exist in the encounter region. Thus the 
contribution from the terms depicted in Figure 1 can be calculated as 

f>T-2te 



t-cnc 



2 t(:»rirl2 
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The HOdA sum rule [IB] changes the sum ^^\Fa\^Ta into / dT. Hence, Eq.(^ is 
rewritten as 

1 r poo /-T— 2tcnc 1 

^ / dsdu / dT / rft^e*W^e-''*/^^e*(^^-^^)^/^«. (33) 

4 i -'2tcnc "'0 "^enc 

To clarify the contributions from the encounter and links, we introduce new 
variables as 



T = ti+t2 + 2t 

enc; 

t =t2. 

Then the expression is simplified as 



(34) 
(35) 



i(eA-eB)t2/TH-bt2/T„ 



xTll ds \ du-^^ 



(36) 



Now we can identify the contributions from the encounter and links: in the above 
expression, the integrals over ti and ^2 are interpreted as the link-contributions, and the 
integral over (m, s) is the encounter-contribution. We expand the expression in tenc and 
extract the constant term, which is expected to survive in the semiclassical limit — >■ 0. 
Then we obtain 

1 



2i[eA - eB + lb) 

We repeat similar calculations for the possibilities listed in Table 1. Summing up the 
results, we find the contribution from the SR pairs 

ryil) _ 2 2 



i{eA- tB + ih) i{eA- tD + ih) 
2 2 



iiec-es + ib) i{ec - eo + ih)' 



(38) 




(b) 




Figure 2. Examples of the parallel Sieber-Richter (aSR) pairs. 



In the evaluation of the second order term, we need to consider the other types 
of pseudo-orbit pairs, two examples of which are shown in Figure 2. These are called 
"parallel Sieber-Richter" (aSR) pairs, following Ref. [H]. In Figure 2(a), an index "A" 
or "C" is assigned to the original orbit (solid curve), and "B" and "D" are to the 



Semiclassical Theory for Universality in Quantum Chaos with Symmetry Crossover 11 



components of the partner pseudo orbit. In Figure 2(6), on the other hand, "A" and 
"C" are assigned to the components of the original pseudo orbit, while "B" or "D" is to 
the partner orbit. Let us calculate the contribution from Figure 2(a). In both of the two 
links, the original orbit and the partner pseudo orbit go in the opposite directions. In 
the encounter, the stretches of the original orbit and those of the partner pseudo orbit 
also go in the opposite directions. Hence, in this case there are constructive interferences 
of the gauge potential in all parts of the pair. The link contribution is given by 



and the encounter yields 



1 POO 1 POO 

Jo Th Jo 



(39) 



J J ^^enc 

In general, when the original (partner) pseudo orbit passes the encounter n+ and 
n_ times (n^ and n'_ times) in the positive and negative directions, respectively, an 
exponential factor exp(mg^j. h t^^clTu) with 

Uenc = ^ |n+ - ?2_ - + n'_ \ (41) 

appears in the encounter contribution [10]. This is a consequence of the stochastic 
behavior of the gauge potential ([9]). Extracting the term surviving in the semiclassical 
limit /i — )■ cxD, we obtain the contribution from Figure 2(a) as 

7^ n ^f_,Ji^A-eB)ti/T„-bh/T„ 

X I— r dt2e'^'^-'^^'^'^"-''^'^^ 
Th Jo 



xT^ [ du [ ds-^e''''/^e'^^^'^-'^--'^^/^"+'^'^'-^^ 

J J \ 'lenc 



2eA - es - eD + i^h 



(42) 



4 (eA - eB + ib){eA - cd + ib) 
A calculation for the diagram in Figure 2(6) is similar and yields 

7(1) -1/^ n df^J{eA-eB)ti/TH-bh/T„\ 

X r dt2e'^^'^-'^)t2/Tn-bt2/TH 

[Th Jo 

xT^ [ du [ ds^LeW'»ei{(eA+^c-eB)/TH+i46}te„c 

- eA + - 6b + i'ib 

4 {eA - eB + ib){ec - es + ib) 
In order to obtain the complete set of the aSR pairs, we need to change the 
assignment of the indices A,B,C and D as well as the directions of the pseudo-orbit 
components. There are 64 diagrams for each of the types in Figure 2(a) and Figure 
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2(b). Thus totally 128 diagrams have to be considered. Summing up the results by 
using " Mathematica" , we find the total contribution from the aSR pairs 

^aSR = 26(eA-ec)(eB-ez3) 

1 



X 



(ec - CD + iWi^A - CD + ih){ec -63 + ib) 
1 



f • (44) 

(e^ -eB + ihY{eA - + ih){ec - + ih) 
Note that, in both of the limits 6^0 and h 00, this expression vanishes: It is finite 
only in the crossover domain. 

3.3. Higher Order Terms 

Let us consider the higher order terms in the periodic-orbit expansion. Suppose that 
the pseudo-orbit pair includes V encounters, as schematically drawn in Figure 3. In the 
figure, the ports (the end points of the encounters) are depicted as black dots. A pseudo 
orbit is formed by connecting the ports by stretches within the encounters and by links 
out of them. In order to see the general structures of the higher order terms, we need 
to count the number of pseudo-orbit pairs, carefully avoiding over countings. 

As we have seen in the calculation of the SR and aSR pairs, we can identify the 
contribution from links and that from the encounters. That is, a link gives a contribution 

' ^- V (45) 

where ^ = if the link is traversed with the same sense in the original and partner 
pseudo orbits, and ^ = 1 if the senses of traversal are opposite. On the other hand, the 
encounter has the contribution 

i{mAeA - rriBeB + mcec - mDen) - bnl^^, (46) 

where m^, ■ ■ ■ , mjj are the number of the stretches belonging to A, B, C, D, respectively. 

Now we estimate the number of over-countings, when there arc V encounters. 
Obviously, the ordering of the encounters is arbitrary. Therefore there are V\ 
overcountings from the ordering of the encounters. Suppose that the jth encounter 
has ij stretches of the original pseudo orbit. Inside the jth encounter, we also overcount 
£j times, since it is arbitrary which port is chosen as the first. Since the stretches are 
no longer required to point from left to right, we can relabel the left side as the right 
side and vice versa. This requires one more overcounting factor 2^. 

Thus the periodic-orbit expansion of is formally written as 



7(1) 



00 ^ 

E E (-1)"^*"° 

V=\ ^i,£2,-"!^i/ pseudo— orbit pairs 

nj=i (2I-) {i(jnA,eA - niB^tB + mc^ec - mD^eo) - bnl^^j} 

niinks(-0{(eAorC " CSori)) + ^/^&} 

(47) 
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where the second sum is over all possible combinations of the number of stretches, 
and the third sum means counting over all the possible pseudo-orbit pairs when 
(£1,4, ■ ■ ■ are given. 




ii = 2 £2 = 4. ' ' ' iv^2 

Figure 3. A schematic picture of a pseudo orbit pair with V encounters. The hnks 
are shown only for the first encounter. The sohd hues depict the original pseudo orbit 
and the dashed lines depict the partner. 



4. Nonlinear Sigma Model in GOE-GUE Crossover Domain 

Let us examine the random matrix theory prediction for the generating function and 
derive an expansion which can directly be compared with the semiclassical formula. For 
that purpose, using the Bosonic variant of the replica trick we analyze the nonlinear 
sigma (NLS) model. Compared with the Fermionic replica and the supersymmetric 
version of the sigma model [19] the Bosonic replica has the disadvantage that it yields 
only the part Z^^^ of the generating function. However we choose it because it is 
technically much simpler, and, anyhow, the complement Z^'^^ can be obtained from 
Z^^^ by simply swapping arguments as ec ^ ez). As shown in the following sections, 
the ensuing perturbative expansion of the NLS model is equivalent to the semiclassical 
periodic-orbit expansion for Z^^\ 

As derived in Appendix A, the NLS model for the generating function is written as 

Z = hm / d[B] exp{5(5)}, (48) 

1 

S{B) = -{eA -eB + ec - en) 

+ '-TiiiAcBB^) - '-TiiisoB^B) - 2sTt{BB^) + 2sTT{Br^B^T3) 

00 00 
+ ^ E [iAciBB^r] - Tr {^^^(Eti?)™}] - 2s ^ Tr(i?i?t)- 

m=2 m=2 
00 00 

-sJ2Yl Tr [{BB^rrsiBB^r'r,] 

m=l m'=l 
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oo oo 

m=l m'=l 

oo 

+ 2s ^ Tr|(55^)™5r3(fi^5)"'fiV3} , (49) 

(m,m')/{0,0) 

where the measure in the integral is defined as d[B] = Yli k^^^BikdlmBik. Iac-, ^bd 
and Ts are 2r x 2r matrices given by 

(-AC = I n^'^ ^ I, ^Ac = diag(eA, ec, ■ - , ed), (50) 

Y U eAc J 

^BD= f ^ |> ^BD = diag(e£,, e^, ■ - , e^;), (51) 

\ U eBD J 

r r 

rs = diag(lX^, -1, -1,^- ■ ■ , -!)• (52) 

The integrations are performed over all the matrix elements of B with the range 
(—00,00). The variable r is the number of replica spaces, and we finally take the 
limit r — )■ 0. 

4.I. Diagonal Term 

Let us calculate the diagonal term ^diag of the generating function. The diagonal 
contribution comes from the quadratic terms of S{B) as 

^diag = e*(^-^-^s-^^+^-)/2 hm / d[B] exp { '-TT{~eAcBB^) - '-TiiiBoB^B) 

-2sTt{BB^) + 2sTr(5r35r3)} . (53) 
In order to explicitly calculate this integral, we decompose the B matrix as 

/ 5++ B+- \ 



where each matrix i?"'^ is an r x r matrix. We note that the symmetry (]A.33ll requires 
the relations 

(i?++)t = - {B-f, (55) 
{B-+y = - {B+-f. (56) 
Using these relations, we calculate fl53|) as 

7,. — ^«(eA-es-ec+':D)/2 

E E { ^(^^^-^ - ^^^^'^) - 2^ + I^M I' 



gi(eA-es-ec+eD)/2 

r-1 

— TT — TT I — TT —71 

X lim 



r-->o 2(eA - esj - 8s HeA-efij [i[<^A - <^d) - as H^A - 
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X 



1 /■ \ (r— l)(r 

-TT — TT I — TT — TT ' 



iiec - es) - 8s i{ec - es) ) { i{ec - eo) - 8s i{ec - cd 



Taking the limit r — )■ 0, we obtain 

7,. — p'i{<^A-<iB-<iC+<^D)/'^ 

(cyi - CD + i8s){ec - + i8s) ) f (e^ - eD){ec - ^b) 



{^A -tB + i8s){tc -tD + i8s) J y{eA- eB){ec - e^) 
This is identical to the semiclassical formula, if 8s is replaced by b. 



■ (58) 



4-2. Expansion of the Generating Function 

The generating function Z is formally expanded as 

^ = 2diag(l + ^off), (59) 



Zoff = lim(^exp i^Te 

I e ' i,e' 



(60) 



where J-> is written as 



= ^Tr{e^c(fi5^)'} - '-Tr{~eBD{B^BY} - 2si:i{BB^Y 

e-i e-i 
- s Tt {{BBY'raiBB^y-'^Ts} - s ^ Tr {(5^5)™r3(5^fi)^"™r3} 

m=l m=l 

+ 2s ^ TT{{BB^rBT3{B^BY-'^~^B^Ts}. (61) 

m=0 

The symbol (...) in Eq.f l60p means an average computed over the Gaussian weight in ( 153|) 
with a normalization factor. Perturbative expansions will be derived by means of the 
Wick theorem on this Gaussian average. Thus the Gaussian average plays a crucial role 
in establishing the relationship between the periodic-orbit theory and random matrices. 
The Gaussian averages of the two-point correlation functions of B are evaluated as 

(BfeY Bi^; ) = - ^M' V'^^.^'W (g2) 

' '''^ / «(eAc,fc -CBD,^) -4s|a-/3|' 

where a' and (3' denote the opposites of the signs a' and /3', respectively, and k, k', i and 
f take the values 1, ■ ■ ■ , r. 

In order to obtain a simple expression for J-"^ in terms of the matrix elements of B, 
we make use of the cyclic property of a trace, i.e., Tr(CiC2 ■ ■ • Ol) = Tr(C2C3 ■ ■ ■ Ci) = 
••• = Tt{OlOi ■ ■ ■ Ol^i). Due to this cyclic property, £ equivalent but seemingly 
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different expressions for J-£ are obtained. Taking the average of tliese expressions, we 
find tlie following compact formula 

~ £ ^kiM y^kaM I ■ ■ ' ^k2e-^_Mt-l^^klMl ) 

ai,---,a2i ki,---,k2t 

^^{^ACM,-x - ^BD,k2,) ~ 2 - a2j-l) > ■ (64) 

_j=l I j=l J 

As clarified below, the squared factor corresponds to the factor n^^^ in the semiclassical 
expression |T0]. In order to evaluate the expansion fj60|) . general correlation functions 
of B are necessary. We are able to calculate them by using the Wick theorem. In the 
next subsection, we obtain the second order term which yields the leading off-diagonal 
contribution. 



1 

X - 
2 



4-2.1. Second Order Contribution Second order contribution in Zos results from the 
quartic terms in B. The quartic terms are simply given by 



,a4=± ki,---,k4=l 



^i[^AC,ki + ^AC,k3 — {^BD,k2 + ^BDm)] " 8-5 ^ 



«! + 0:3 0^2 + ^4 



X 



We employ the Wick theorem to write this in terms of the two-point correlation functions 
( l62l) and ( l63l) . That is, a four-point correlation function is expanded as 

j^ai,a2 



■'ki,k2 y k3,k2 ) k3,ki 



,04 \* 

fe4 ) 



+ 
+ 



71)01,02 003,04 



pOi,a2 / p03,02\* 



k\,k2 



pOl,02 

^kxM 



I 7503,02 \* 

yBk^M ) 



7^03,04 

fc3,fc4 



^BZ 



,04 \* 

k^ ) 



( 00l,04\ * 



( pOi,04\* 



003,02 / 7003,04 \* 

^3,^2 k3,k4, ) 



In the periodic-orbit theory, the SR and aSR pairs give the corresponding contribution 
to these terms. Let us consider the contribution from the first term 

5Z 5Z (^) ^(2eAc,A,-i -'iesDM) 

01,02 fel,fc2 

1 If 1 



K^ACM ~ ^bdm) - 4s|ai - j yi{eAc,ki - ^bdm) - 4s|ai + 02 

(66) 

This expression have 16 terms after taking the replica limit r — t- 0, which is consistent 
with the the number of the SR pairs in the semiclassical theory. As easily checked, this 
reproduces the result in (!38|) . if 8s is replaced by h. The second term generates the 
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contribution corresponding to that of the aSR pairs of the type in Figure 2(a): 



i{2eAc,ki - i^BDM + ^bdm)} -8s [ai 
1 



a2 + «4 



X 



ai—a2 
2 I 



a\ — Q4 
I 2 I 



(67) 

This has 64 terms after taking r — )■ 0. Finally the third term gives the contribution 
corresponding to that of the aSR pairs of the type in Figure 2(6): 

01,03,02 ki,k3,k2 

ne^cfci + CACfcg - 2eBD,fc2) - 8s I 02 I 



^(cyicfci — e_BD,fc2) ~ 8s "^^"^ J 1 "i-i^ACks — ^BD,k2) " 8s 



03—02 
2 



(68) 



Eqs. (!67|) and ( |68|) totally yield 128 terms, and reproduce the result in (l44l) . 



z^.,?. Relation between Semiclassical Formulas and Random Matrix Results 

Let us consider how the semiclassical formulas are related to random matrix results. In 
the expansion ( 160|) . the factor J^^ corresponds to the ^-encounter, while the l^-th order 
terms in J^i are related to the terms with V encounters. The factors B and B* represent 
the ports of the encounters. The contraction lines in the Wick theorem correspond to the 
links. Let us consider the factors -B^"^^^^ and {B^^'^^'^^^)* in J^i and the corresponding 
ports. Note that a2i+i and k2i+i are identified with ai and ki, respectively. The ports 
are connected inside the ^-encounter as follows. If i + 1 = j, the original pseudo orbit 
connects the ports. On the other hand, if i = j ' + 1, the partner connects them. The 
superindex a on the left (right) side specifies the direction of the partner (original) 
pseudo orbit. The subindex k on the left (right) side determines the assignment of the 
index B oi D {A or C). Using these rules, we can easily interpret Eq.(l66l)-(l68l) as the 
SR and aSR pairs. 

As an example, let us examine the following case 




where the lines imply contractions. This example consists of two encounters with 2 
and 3 stretches. In Figure 4, the solid curves depict the original pseudo orbit, and the 
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dashed curves depict the partner. The port B is located on the left side and B* is on 
the right side of each encounter. In the figure, only the port indices and the connections 
among them are written. The thin curves (lines) express the connections inside the 
encounters, and the bold curves express the links generated by taking the contractions. 
Inside each encounter, the ports are connected by following the explained rules. The 
directions of the links and the encounter stretches are also shown in the figure . If the 
link connects B and B*, then their left superindices (e.g., ai and a'^) are the same and 
so are their right superindices (e.g., and a'^). On the other hand, if the link connects 
B and B, or B* and B*, their left (right) superindices become the opposite. These 
properties originate from ( l62l) and ( l63l) . The signs a can consistently be interpreted as 



the directions of the links. In Figure 4, the signs a are ai = 



+, 



04 



an 



and a2 = +■ 



In the semiclassical theory, for each topological structure of the pseudo-orbit pairs, 
we have to consider all the possible ways of assigning either A or C to the components 
the original pseudo orbit, and assigning either B or D to those of the partner. We see 
that each component accompanies a factor —1, if C or D is assigned to it. In the replica 
treatment of random matrices, on the other hand, if C or D is assigned, the summation 
over all the possible choices of the indices leads to a factor r — 1, which goes to —1 in 
the limit r — > 0. Hence the sign finally becomes the same as that in the semiclassical 
formula. 



_T__ 

Oil 

ki 
- 1" 
„ i_ 



Q^2 
k2 



"1" 



0^2 
k2 



I 

I 



«2 
^2 



«3 



^3 



k 



k/i 



cy.1 
ki\ 



Qi4 

k4 



_L _ 



«4 



etc 



0^4 



5 1 



«6 
^6 



i 



«6 

kc 



Figure 4. An example which consists of two encounters. Solid curves depict 
the original pseudo orbit, while the dashed curves depict the partner. The links 
correspond to the contractions in the Wick theorem. The directions of the links and 
the encounter stretches are also shown. In this case, ai = = a'^ — a'l — —Ci-^ — +, 
a4 = — ag = —o!a = — , o!o = — , and a2 = +■ 
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Finally, we discuss how many times the same structure appears, which gives the 
overcounting factor. Let us consider the ^-encounter term Eq. flMj) . and suppose that 
the contractions are already given. 

First we change the indices as 

), (69) 
) (70) 

for m = For each value of p = 1, ■ ■ ■ , £, we find the same pseudo orbit pairs. 

This is related to the fact that we can arbitrarily choose one of the i ports as the first. 
Next we exchange the left and right ports in addition to the transformation 

(«!, ^i) ("3, h), (71) 

(as, ks) (ai, fci), (72) 

{am+3,km+3) {a2e-i+m, he+i-m) , m > 1, (73) 

{a2e+i-m, k2e+i-m) -> (as+m, ^s+m), m>l. (74) 

That is, changing the layout in Figure 4, we put the indices of B and B* on the right 
and left sides, respectively. Then we recover the original pair of the pseudo orbits. This 
procedure yields two identical pairs. This related to the fact that we can relabel the left 
side of an encounter as the right side and vice versa in the semiclassical argument. 

Thus we found that the overcounting factor was in total 2i for each structure of 
the pseudo-orbit pairs. 

4.4- Diagrammatic Expansion of the Generating Function 

Now we are able to write down a diagrammatic expansion of Zos- Using (l60ll and the 
Wick theorem, we obtain a formal expression 



Z, 



off 



X 



V=l ' lll2,-/v 
oo ^ 

V=l ^i/2,---/v contractions k,a 

X {l/2ij){i{mA^eA - me^eB + mc^ec - mD^eo) - 8s nl^^ j} 
jJl (-l){^(eAorC - efiorl?) - 8SyU} 



(75) 

where the symbol ^ ' stands for the sum over all allowed k and a after taking the 
contractions. The corresponding values of niAj, ■ ■ ■ , moj, nencj, A* and the assignments 
of e^orC) CfiorD are diagrammatically determined for each term. Note that a factor 
is included in order to make up for the overcountings discussed above. We conjecture 
that this diagrammatic expansion is identified with the periodic-orbit expansion (H7!) 
derived in 53. 
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5. Summary 

For a classically chaotic system in a weak magnetic field, the generating function of 
the spectral correlation was studied by using the periodic-orbit theory. The diagonal 
and leading off-diagonal terms in the periodic-orbit expansion were explicitly calculated. 
Both the non-oscillatory and oscillatory terms were in agreement with the prediction of 
the random matrix theory. Moreover we clarified the general structures of the higher 
order terms. In order to see the correspondence to the result of the random matrix 
theory, we analyzed the nonlinear sigma (NLS) model in the crossover domain between 
the GOE and GUE classes. It was conjectured that the diagrammatic expansion of the 
NLS model could be identified with the periodic-orbit expansion. 



Acknowledgements 

The authors are grateful to Prof. Fritz Haake and Dr. Stefan Heusler for stimulating 
discussion and their continuous interest. They also thank Prof. Martin Zirnbauer for 
providing his original result [15] before publication. This work is partially supported 
by Japan Society for the Promotion of Science (KAKENHI 20540372) and by the 
Sonderforschungsbereich TR 12 of the Deutsche Forschungsgemeinschaft. 



Appendix A. Random Matrix Theory in GOE-GUE Crossover Domain 
Appendix A.l. Derivation of the Nonlinear Sigma Model 

In this Appendix, we derive the nonlinear sigma model in the crossover domain between 
the GOE and GUE classes of random matrices. The probability distribution of the 
random matrix H reads 

N 



P{H)dH = Xexp 



N 



4{((To -fTu)e-2- + au} 



j<i 



N 



0-2 



dH, 



(A.l) 



2(Tu (1 - e 

where A/" is a normalization factor. Throughout this Appendix, we use the symbol 
A/" for normalization factors, even if their values are different. The distribution (lA.ip 
reproduces the GOE distribution Pgoe in the limit r — > 0, and the GUE distribution 
Pgue in the limit r — )■ oo, respectively: 

P{H) Pgoe(^) 

TV 



oc exp 



A^ 



4ao 

P{H) ^ PGUE(i^) 

N 



r2. 



oc exp 



N 



J j<e 



r ^ 



4(Tt 



r — !■ oo 
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We use the Bosonic replica trick to calculate the generating function defined by (1) and 
derive the nonlinear sigma model (NLS). Here we use the notation Z to discriminate 
the NLS expression from the random matrix expression (3). The generating function is 
formally written in terms of the vectors V'^, V^b, il^c,, and il^Ot, {k = 1, • • • , r — 1), each of 
which consists of elements: 

Z = lim y d[V^]e*'^^^-^'^ 

/ / r-l r-1 

X ( cxp I -iij\HlpA - ?■ Yl ^Ck^^Ck + ii^B^^B + i Y '^Dk^^Dk 
\ \ k=l k=l / 

(A.2) 

where ip is a vector composed of ipA, ipc^ {k = 1, ■ ■ ■ , r — 1), t/'s, and tpn,^ (fc = 1, ■ • • , r — 1). 
The symbol (...) denotes the random matrix average over the measure (A.l), and E and 
L are 2rN x 2rN matrices 

7 1 1 1 



E = diag(£;+l, E+1, E+1, E-1, E'l, E-1) , (A.3) 



r 



L = diag(l, 1, ■ ■ ■ , 1, -1, - . (A.4) 

Here 1 is the N x N identity matrix. A random matrix average (...) can be easily 
evaluated, and we obtain 

Z = lim| d[^]e^^^^^^exp (-^^Va^ + ^^TVft^^ , (A.5) 

where = (ctq — a[/)e^^'^ + and /x^ = au{l — e"^"^). The matrices a and b are N x N 
matrices, whose elements are given by 

ae,m = Yl %A'''^n,m + Y , (A.6) 

r) T) 

Here i] is one of the indices ^4, Ci, • • • , C^-i, B,Di, - ■ ■ , -Dr-i- Let us define hJ'^'^' as the 
(r?, T]') element of the diagonal matrix 



A = diag(LL---,l,-l,-l,---,-l). (A.8) 
Next we introduce the vector and ^ which have 4r components as 

Then Tra^ and Tr6^ can be written in terms of these vectors as 

TVa^ = 4 ^ ^iM'^rJfmM'^i, (A.IO) 



Trft^ = - 4 J] ^,M-^J>^M'^,, (A.ll) 
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where M and M are defined as 

By introducing 4r x 4r matrices A and A as 

N N 
m=l m=l 

we can write the generating function in the form 

J rf[^]e«'^^^^'^e-^(^^'+^'^^'), (A.14) 

where u'^ — (fx/X)^. Here and hereafter we omit the symbol hmr_>o- In order to write A 
and A in a unified way, we introduce A as 

A = b+A + fe-fg^fg, (A.15) 

6± = ^ (\/rT^ ± Vl^^ , (A.16) 

where fg is 

2r 2r 

fs = diag(17v^,-l,-V--,-l). (A. 17) 

Notice that A = Af^, and Tr^^ = TrA'^ + z/^Tr[(74f3)^]. Thus, the generating function 
is written in terms of the new matrix A as 



We then employ the Hubbard- Stratonovich transformation to get rid of the 
term quadratic in A and consequently quartic in ■0. Using the relation 
/ exp (^TrQ^) = 1, where Q is anti-Hermitian, we can readily derive the 
expression 

j d[Q\ j d[^l^\ exp + ^Tr<5^ - ATr(^Q) | . (A.19) 

Noting that the term 1i{AQ) is written as Tx[AQ) = TY{A{b+Q + b-f^Qfs)}. Hence, 
defining new matrix Q as 

Q = b+Q + b.fsQfs, (A.20) 

we obtain 



Af J d[Q] J ci[^] exp 1^ J2 ^rnM (^F + iAg) + jTiQ^j 



(A.21) 

Here F is a 4r x 4r matrix defined as 

r— 1 r—1 

F — diag(£^A, Ec, • • • , Eq, Eb, Ed, • • • , Ed, 
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In order to obtain Eg. ([OB . the relation i/j^fLEip = ^ip^ LEip + ijj^ ELip) = 
\ [ip'^LEip + iIFLEiI)*) = J2m=i ^m.MF'^rn was used. Evaluating the ip-integral, we 
reduce Eg- lTOT]) to 



Z=M j d[Q] exp I jTVQ^ - y Tr ln(F + i\Q) 

We set the average energy E to the origin (the band center), and write F as 
1 



(A.23) 



(A.24) 



27rp 

1 — 1 r— 1 r— 1 r— 1 

e = diag(e^, ec, • • • , ec, cb, Ci?, ■ ■ ■ , e^, ^a-, ^c, ■ ■ ■ , ^c, ^b, ^d, ■ ■ ■ , ^d) , 

(A.25) 

where the average density is p = N/IttX). Now we need to consider the effect of u. We 
treat it as a perturbation [20] and expand 6+ and 6_ with respect to as 6+ ~ 1 and 
6_ ~ Setting u = 0{N~^^'^), we find a dominant contribution to the generating 
function 



Z=Af I d[Q] exp { -yTrln ( 1 + 



e + 



iXQ 



X exp 



iV<;iTrQ2_lTrln(zAg) 



(A.26) 



The next step is the saddle point analysis in the limit N oo. We find the saddle 
point equation 



5 



^■Trg2 _ iTrln(zAg) }> = . 



This equation yields a saddle point satisfying 



Ql 



1 . 



(A.27) 



(A.28) 



As in the case of the GOE, the saddle point satisfies the symmetry relation 

= . (A.29) 



Here cxi is the Pauli matrix 



0, 1 

1, 

we obtain the nonlinear sigma model 



. Expanding the logarithmic part of Eq. ( 1A.261) . 



Z=M j d[Q,] expj-yTrln ^1 + 



= A/'y rf[Q*]exp 

= j c?[Q*]exp 
Here we introduced a parameter s = Nv"^. 



i Nv"^ 
-Tr(eQ.)-^Tr{(af3r} 



^Tr(6Q.)-^Tr{(Q.f3f} 



(A.30) 
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Appendix A. 2. Parametrization 

We further consider the exponential factor in Eq. ( 1A.30I) . Let us denote the exponent 
in Eg. (1X301) by 

5(g.) = ^Tr(eg.) - ^Tr {{Q^f,)'} , (A.31) 
As in the GOE case, we parametrize by using the matrices B and as 



0, B \ \ 1, \ I , / 0, B 

fit, j 0, -1 j r 



n 

0, B \ \ I 1, \ / 0, B 



1, \^ 

0, -1 / ^ 

' / n=0 



(A.32) 



In this representation, the symmetry relation of Qs (]A.29|) is satisfied by requiring 

= -aifiV- (A.33) 

From a straightforward calculation, we obtain the expression of S{B) ( = ) 
in the form 

S{B) = -{ea -eB-ec + e^) 

+ '-TiiiAcBB^) - '-TiiiBoB^B) - 2sTt{BB^) + 2sTt{BtsB^ts) 

oo oo 

+ ^ E {iAciBBY'} - Tr {eBD(fi^fi)™}] " 2s Tr(i?i?t)- 

m=2 m=2 
oo oo 

-sJ2Yl Tr{(i?fit)-r3(fii?t)-'r3} 

m=l m'=l 

oo oo 

-sJ2J2 Tr{iB^Br%{B^By^'r,] 

m=l m'=l 

oo 

+ 2s J2 T^{{BBY'BT3{B^B)"''Bh3^ , (A.34) 

m,m'=0;(m,m')7^{0,0) 

where e^csD and are 2r x 2r matrices defined in ( 15U]) . and 
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